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We discuss several scenarios for the creation of nonlocal spin-entangled electrons which provide 
a source of electronic Einstein-Podolsky-Rosen (EPR) pairs. Such EPR pairs can be used to test 
nonlocality of electrons in solid state systems, and they form the basic resources for quantum infor- 
mation processing. The central idea is to exploit the spin correlations naturally present in super- 
conductors in form of Cooper pairs possessing spin-singlet wavefunctions. We show that nonlocal 
spin-entanglement in form of an effective Heisenberg spin interaction is induced between electron 
spins residing on two quantum dots with no direct coupling between them but each of them be- 
ing tunnel-coupled to the same superconductor. We then discuss a nonequilibrium setup with an 
applied bias where mobile and nonlocal spin-entanglement can be created by coherent injection of 
two electrons in a pair (Andreev) tunneling process into two spatially separated quantum dots and 
subsequently into two Fermi-liquid leads. The current for injecting two spin-entangled electrons into 
different leads shows a resonance and allows the injection of electrons at the same orbital energy, 
which is a crucial requirement for the detection of spin-entanglement via the current noise. On the 
other hand, tunneling via the same dot into the same lead is suppressed by the Coulomb blockade 
effect of the quantum dots. We discuss Aharonov-Bohm oscillations in the current and show that 
they contain h/e and h/2e periods, which provides an experimental means to test the nonlocality 
of the entangled pair. Finally we discuss a structure consisting of a superconductor weakly coupled 
to two separate one-dimensional leads with Luttinger liquid properties. We show that strong corre- 
lations again suppress the coherent subsequent tunneling of two electrons into the same lead, thus 
generating again nonlocal spin-entangled electrons in the Luttinger liquid leads. 



I. INTRODUCTION 

The electron has charge as well as spin. While the 
control of charge is well established in electronic systems 
and devices, the spin combined with transport only re- 
cently attracted interest, both from a fundamental point 
of view and also for applications in electronicaiEl. The 
idea of using the spin and its coherence for electronic de- 
vices has received strong experimental suppprtuLl show- 
ing unusually long spin dephasing timea3ij in semi- 
conductors (approaching microseconds), the injection of 
spin-polarized currents from a magnetic to a nonmag- 
netic semiconductorBQ, as well as phase-coherent spin 
transport over distances of up to 100 /mua til. Besides 
the improvement of conventional devicesM due to spin, 
e.g. in magnetic read out heads, non-volatile memories 
etc., electron spins in quantum confined nanostructures 
like semiconductor quantum dots, atoms, or molecules 
can be usecLas a quantum bit (qubit)B fQx._q.uantum 
computingmia and quantum communicatiorMEJ, where 
a radically new design of the necessary hardware is re- 
quired. In this article we focus on the creation of spin- 
entangled electrons-electronic Einstein-Podolsky-Rosen 
(EPR) pairalj-exploiting the spin-correlations naturally 
present in s-wave superconductors, where the electrons 
are correlated as Cooper pairs with spin-singlet wave- 
functions. We note that such EPR pairs represent the ba- 
sic resources for quantum communication!!!! schemes like 



dense coding, teleportation, or, more fundamentally, for 
testing nonlocality via Bell inequalitiesLj, which would be 
particularly interesting to implement for massive paijti.- 
cles such as electrons in a condensed matter systemE£lli3. 

We describe in the following a few proposals for cre- 
ating nonlocal spin-entanglcmcnt between two electrons. 
We consider first two quantum dots which are coupled to 
the same superconductors (but not among themselves). 
In equilibrium, this coupling .then induces an effective 
tunable Heisenberg interaction!^ among the two electrons 
on the_dots which could be used to implement two-qubit 
gateaij. We then consider a nonequilibrium situation 
where now Cooper pairs can be transported by means 
of an Andreev (pair tunneling) process from the super- 
conductor via two quantum dots into different normal 
Fermi liquid leads by appkung a bias between the su- 
perconductor and the leadsEH. The quantum dots in the 
Coulomb blockade regime are needed to mediate the nec- 
essary interaction between the two tunneling electrons 
so that they preferably tunnel to different dots and then 
subsequently to different outgoing leads, thereby main- 
taining their spin-singlet state. Such a setup then works 
as an entangler for electron spins, satisfying all necessary 
requirements, within a parameter regime of experimen- 
tal accessebility, to detect-the spin-entanglement in the 
leads via the current noisetil. We finally discuss a further 
realization for an entangler where the necessary interac- 
tion to separate the two electrons is provided by strong 
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correlations in one-dimensional leads (such as nanotubes) 
with Luttinger liquid properties. 

We refer to related worktS which makes also use of 
Andreev tunneling but with a transparent superconduc- 
tor/normal interface. The electrons move from the su- 
perconductor into a normal fork-shaped wire without 
Coulomb blockade behavior. The electrons are separated 
via energy filters so that the electrons enter their cor- 
responding leads at different orbital energies. However, 
due to the transparent interface the partners of different 
Cooper pairs are not separated in time and space which 
is needed to identify the spin-entangled partners. Subse- 
quent work-ip-similar direction uses multiterminal hybrid 
structureal3'E2l. 



II. CREATION AND DETECTION 
OF NONLOCAL 
SPIN-ENTANGLEMENT IN A DOUBLE DOT 

We consider a parallel double quantum dot structure. 
The quantum dots contain one electron spin each and 
we assume that there is no direct coupling between the 
dots but each dot is tunnel- coupled to the same s-wave 
superconductors with a tunnel amplitude T (see Fig. 1) 
on both sides of the dots. The two superconductors are 
held at the same chemical potential. The s-wave super- 
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FIG. 1. Upper panel: Sketch of the superconductor-double 
quantum dot-superconductor (S-DD-S) nanostructure. There 
is no direct coupling between the dots. Lower panel: 
Schematic representation of the quasiparticle energy spectrum 
in the superconductors and the single-electron levels of the 
two quantum dots. 

conductor favors an entangled singlet-state on the dots 
(like in a Cooper pair) and further provides a mechanism 
for detecting the spin state via the Josephson current. It 
turns out that in leading order oc T 4 the spin coupling 



is described by a Heisenberg HamiltoniarJl3 ! EJ with an 
antiferromagnetic exchange coupling J 

ff cff « J(l + cos^) (s a -S b -^j , (2.1) 

where J « 2T 2 /e, and the energy of the dot is e below the 
chemical potentials of the superconductors. Here, (p is 
the average phase difference across the superconductor- 
double-dot-superconductor (S-DD-S) junction. We can 
modify the exchange coupling between the spins by tun- 
ing the external control parameters T and ip so the device 
can act as a two-qubit quantum gate needed for quantum 
computing where the-,electron spin is the fundamental 
computational unitErEa. Furthermore, the entangled spin 
state on the dot can be probed if the superconducting 
leads are joined with one additional (ordinary) Joseph- 
son junction with coupling J' and phase difference 8 into 
a SQUIEkring. The supercurrent Is through this ring is 
given byEj 

, _ J sin(0-27r/) + (J7</)sin0, singlet, , . 
s/ J ~ \ (J'/ J) sin 9, triplets, 1 > 

where Ij — 2eJ/Ti and / = $/$o with <E> being the mag- 
netic flux threading the SQUID-ring and $ = Tic/2e is 
the flux quantum. Measurement of the spin- and flux- 
dependent critical current I c = maxg{|/s|} probes the 
spin state of the double dot. This is realized by biasing 
the system with a dc current / until a finite voltage V 
appears for |/| > ljt5- 

III. ANDREEV ENTANGLER 

We now consider a nonequilibrium situation where the 
superconductor is tunnel-coupled to two quantum dots in 
the Coulomb blockade regime which are further coupled 
to two normal Fermi liquid leadalU, see Fig. 2. Applying 
a bias between the superconductor and the leads than re- 
sults in a stationary current of spin-entangled electrons 
via Andreev tunneling and resonant transport from the 
superconductor to the leads. The quantum dots are used 
to mediate the interaction necessary to separate the two 
spin-entangled electrons originating from a Cooper pair. 
The amount of spin-entanglement in the outgoing cur- 
rent could be tested via noise measurementsE-i We have 
shown that current-current correlations (noise) are en- 
hanced if the injected electrons are singlets due to bunch- 
ing behavior whereas the noise is suppressed in the case of 
spin triplets due to antibunching behavior. For such noise 
measurements, which are based on two-particle interfer- 
ence effects, it is absolutely crucial that both electrons, 
coming from different leads, possess the same orbital en- 
ergy. 
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FIG. 2. The entangler setup: Two spin-entangled electrons 
forming a Cooper pair can tunnel with amplitude Tsd from 
points ri and r2 of the superconductor, SC, to two dots, Di 
and D2, by means of Andreev tunneling. The dots are tun- 
nel-coupled to noninteracting normal leads Li and L2, with 
tunneling amplitude Tdl- The superconductor and leads are 
kept at chemical potentials [is and /iz, resp. 



IV. QUALITATIVE DESCRIPTION OF THE 
ANDREEV ENTANGLER 

We first provide a qualitative description of the en- 
tangler and its principal mechanism based on Andreev 
processes and Coulomb blockade effects and also specify 
the necessary parameter regime for successful transport 
of the initial spin-entanglement of the Cooper pairs to the 
outgoing leads. In subsequent sections we then introduce 
the Hamiltonian and calculate the stationary current for 
two competing transport channels which is followed by a 
discussion of the results. 

We consider an s-wave superconductor where the elec- 
itpns form Cooper pairs with singlet spin-wavefunctions 
E3. The superconductor, which is held at the chemical 
potential us, is weakly coupled by tunnel barriers to two 
separate quantum dots D\ and D2 which themselves are 
weakly coupled to Fermi liquid leads L\ and L2, resp., 
both held at the same chemical potential ji\ = ji% = fJ-l- 
The corresponding tunneling amplitudes between super- 
conductor and dots, and dot-leads, are denoted by Tsd 
and Tdl, resp. which, for simplicity, we assume to be 
equal for both dots and leads. 

By applying a bias voltage Afi = /4<j — /J-i > transport 
of entangled electrons occurs from the superconductor via 
the dots to the leads. In general, the tunnel-coupling of 
a superconductor to a normal region allows for coher- 
ent transport of twQ electrons of opposite spins due to 
Andreev tuanelingE3, while single-electron tunneling is 
suppressed!^ in the regime A > A/i, ksT, where A is the 
energy gap in the superconductor and T is the temper- 
ature. The gap A is the minimum energy to break up 
a Cooper pair into a quasiparticle in the superconductor 
and an electron in the normal region due to tunneling. 
According to the energy-time uncertainty relation, h/A 
then defines the time delay between the two coherent 



tunneling steps in the Andreev process. In the present 
setup, we envision a situation where the two electrons 
are forced to tunnel coherently into different leads rather 
than both into the same lead. This situation can be 
enforced in the presence of two intermediate quantum 
dots which are assumed to be in the Coulomb blockade 
regimeE-3 so that the state with the two electrons being on 
the same quantum dot is strongly suppressed, and thus 
the electrons will preferably tunnel into separate dots and 
subsequently into separate leads (this will be quantified 
in the following). 

The chemical potentials e\ and 62 of the quantum dots 
can be tuned by external gate voltageaHj such that the 
coherent tunneling of two electrons into different leads is 
at resonance if £1+62 = 2/i£, see Fig. 3. This current 




e 2 



FIG. 3. The energy situation of the superconductor with 
chemical potential /is and the two dots 1,2 with chemical 
potentials ei,£2- Transport of the two members of a Cooper 
pair with energy 2/is from the superconductor to different 
outgoing leads with chemical potential /ii < ei (not drawn) is 
at resonance if ei + £2 = 2/is- 

resonance condition reflects energy conservation in a tun- 
nel process of a Cooper pair with energy 2/tg from the 
superconductor to the dots 1,2 (one electron on each dot) 
with chemical potentials e 1; e^ and requires that the reso- 
nant dot levels have to be adjusted such that one is above 
Us and the other (by the same amount) below [is- This 
is very similar to the more familiar picture of Andreev 
reflection at a superconductor/normal interface. There 
an electron on the normal side of the junction, and with 
energy e above [is > is back reflected as a hole with energy 
e below [is by the simultaneous creation of a Cooper pair 
in the superconductor. In that sense the empty dot level 
below [is can be considered as the hole and the empty dot 
level above [is as the empty electron state. In contrast, 
we will see that the current for the coherent tunneling of 
two electrons via the same dot into the same lead is sup- 
pressed by the on-site Coulomb U repulsion of a quantum 
dot and/or by the superconducting gap A. 

Next, we introduce the relevant parameters describing 
the proposed device and specify their regime of interest. 
First we note that to avoid unwanted correlations with 
electrons already on the quantum dots one could work 
in the cotunneling regimea where the number of elec- 
trons on the dots are fixed and the resonant levels e;, 
I = 1, 2 cannot be occupied. However, we prefer to work 
at the particular resonance ti ~ [is, since then the to- 
tal current and the desired suppression of tunneling into 
the same lead is maximized. Also, the desired injection 
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of the two electrons into different leads but at the same 
orbital energy is then achieved. In the resonant regime, 
we can avoid unwanted correlations between tunneling of 
subsequent Cooper pairs if we require that the dot-lead 
coupling is much stronger than the superconductor-dot 
coupling, i.e. \Tsd\ < \Tdl\, so that electrons which en- 
ter the dots from the superconductor will leave the quan- 
tum dots to the leads much faster than new electrons 
can be provided from the superconductor. In addition, 
a stationary occupation due to the coupling to the leads 
is exponentially small if A/i > kgT, T being the tem- 
perature and ks the Boltzmann constant. Thus in this 
asymmetric barrier case, the resonant dot levels e/ can be 
occupied only during a virtual process (see also Section 

Next, the quantum dots in the ground state are al- 
lowed to contain an arbitrary but even number of elec- 
trons, No = even, with total spin zero (i.e. antiferro- 
magnetic filling of the dots). An odd number Njj must 
be excluded since a simple spin-flip on the quantum dot 
would be possible in the transport process and as a re- 
sult the desired entanglement would be lost. Further, 
we have to make sure that also spin flip processes of the 
following kind are excluded. Consider an electron that 
tunnels from the superconductor into a given dot. Now, 
it is possible in principle (e.g. in a sequential tunneling 
processed) that another electron with the opposite spin 
leaves the dot and tunnels into the lead, and, again, the 
desired entanglement would be lost. However, such spin 
flip processes will be excluded if the energy level spacing 
of the quantum dots, Se, (assumed to be similar for both 
dots) exceeds both, temperature k B T and bias voltage 
A/i. A serious source of entanglement-loss is given by 
electron hole-pair excitations out of the Fermi sea of the 
leads during the resonant tunneling events. Since then a 
simple spin flip on the dot would be possibkjdue to the 
coupling to the leads. However, we showedO that such 
many-particle contributions can be suppressed if the res- 
onance width 7/ = 2itvi\Tdl\ 2 is smaller than A/z (for 
ei ~ fis)i where vi is the density of states (DOS) per 
spin of the leads at the chemical potential /i/ . 

To summarize, the regime of interest where the coher- 
ence of an initially entangled Cooper pair (spin singlet) 
is preserved during the transport to the leads is given by 



A, U, Se > A/x > 7/, k B T, and 7/ > 75 



(4.1) 



As regards possible experimental implementations of 
the proposed setup and its parameter regime, we would 
like to mention that, typically, quantum dots are made 
out of semiconducting heterostructures, which satisfy 
above inequalitiescJ. Furthermore, in recent experi- 
ments, it has been shown that the fabrication of hybrid 
structures with semiconductor and superconductor being 
tunnel-coupled is possibleOcj. Other candidate materi- 
als are e.g. carbon nanotubes which also show Coulomb 
blockade behavior with U and Se being in the regime of 
interest hereEH. The present work might provide further 
motivation to implement the structures proposed here. 



Our goal in the following is to calculate the stationary 
charge current of pairwise spin-entangled electrons for 
two competing transport channels, first for the desired 
transport of two entangled electrons into different leads 
(Ji) and second for the unwanted transport of both elec- 
trons into the same lead (12)- We compare then the two 
competing processes and show how their ratio, I\/ 12, de- 
pends on the various system parameters and how it can 
be made large. An important finding is that when tun- 
neling of two electrons into different leads occurs, the 
current is suppressed due to the fact that tunneling into 
the dots will typically take place from different points ri 
and r 2 on the superconductor (see Fig. 1) due to the spa- 
tial separation of the dots D\ and Z?2- We show that the 
distance of separation Sr = |ri — r2 1 leads to an exponen- 
tial s uppr ession of the current via different dots if Sr > £ 



(see 



3)), where £ = vf/ttA is the coherence length 
In the relevant regime, Sr < £, how- 



of a Cooper pair, 
ever, the suppression is only polynomial in the parameter 
kpSr, with kp being the Fermi wavevector in the super- 
conductor, and depends sensitively on the dimension of 
the superconductor. We find (see Section [x]) that the 
suppression is less dramatic in lower dimensional super- 
conductors where we find asymptotically smoother power 
law suppressions in kpSr. On the other hand, tunneling 
via the same dot implies 5r = 0, but suffers a suppression 
due to U and/or A. The suppression of this current is 
given by the small parameter (^1/U) 2 in the case U < A, 
or by (7;/A) 2 , if U > A as will be derived in the follow- 
ing. Thus, to maximize the efficiency of the entangler, 
we also require kpSr < A/ji, U/~fi. 

Finally, we will discuss the effect of a magnetic flux on 
the entangled current in an Aharonov-Bohm loop, and we 
will see that this current contains both, single- and two- 
particle Aharonov-Bohm periods whose amplitudes have 
different parameter dependences. This allows us to dis- 
tinguish processes where two electrons travel through the 
same arm of the loop from the desired processes where 
two electrons travel through different arms. The rela- 
tive weight of the amplitudes of the two Aharonov-Bohm 
periods are directly accessible by flux-dependent current 
measurements which are then a direct probe of the de- 
sired nonlocality of the entangled electrons. 



V. HAMILTONIAN OF THE ANDREEV 
ENTANGLER 

We use a tunneling Hamiltonian description of the sys- 
tem, H = Ho + Hp, where 

H q = H s +J2Hdi+J2 H ^ 1=1,2. (5.1) 



Here, the sjrperconductor is described by the BCS- 
HamiltonianB H s = J] k(r EkJuoTk*, where a =t,|, 
and the quasiparticle operators 7ko- describe excitations 
out of the BCS-groundstate |0)s defined by 7k<r|0)s = 0. 
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They are related to the electron annihilation and cre- 



ation operators-, Ck 
transformation!^! 



and 



through the Bogoliubov 



CkT = "k7kT + ' y k7-ki 
c-ki = «k7-ki - Wk7kT 



(5.2) 



where u k = (1/V2)(1 + 6c/^k) 1/2 and Uk = (1/^3) t 1 - 
Ck/^k) 1 ^ 2 are the usual BCS coherence factorsca, and 
£k = e k ~ /is is the normal state single-electro n energy 
counted from the Fermi level /is, and = + A 2 
is the quasiparticle energy. We choose energies such that 
/is = 0. Both dots are represented as one localized 
(spin-degenerate) level with energy ei and is modeled by 
an Anderson-type Hamiltonian Hdi = ei J2a A,„di a + 
Uni^nii, I = 1,2. The resonant dot level e/ can be 
tuned by the gate voltage. Other levels of the dots do 
not participate in transport if 5e > Ap > fc^T, where 
Ap = — fii, and pi is the chemical potential of lead 
I = 1,2, and 5e is the single particle energy level spac- 
ing of the dots. The leads I = 1,2 are assumed to be 
nonintcracting (normal) Fermi liquids with Hamiltonian 

Hli = Sko- e k a 3ko- a Zko-- Tunneling from the dot I to the 
lead / or to the point r; in the superconductor is described 
by the tunnel Hamiltonian Hj- = H$d + H dl with 



H sd = TsodiMn) + H.c, 

la 

H D L = ^ T DLal ka di a + H.c. . 



(5.3) 
(5.4) 



Iko 



Here, ip a {vi) annihilates an electron with spin a at site r;, 
and d\ a creates it again (with the same spin) at dot I with 
amplitude T$d- Vv( r z) is related to Ck CT by the Fourier 
transform \j} a {ri) — e lkri c ka . Tunneling from the dot 
to the state k in the lead is described by the tunneling 
amplitude Tul- We assume that the k-dependence of 
Tdl can be safely neglected. 



number which we take to zero at the end of the calcula- 
tion. Finally, pi is the stationary occupation probability 
for the entire system to be in the state \i). The T-matrix 
T(si) can be expanded in powers of the tunnel Hamilto- 
nian Ht, 



1 



^ v ei + ir) - H { 

n—l 



■H T y 



(6.3) 



where the initial energy is Ei — 2 (is = 0. We work 
in the regime defined in Eq. (4.1), i.e. 7/ > 75, and 
A,U,Se > A(i > 7;,fcflT, and around the particular 
resonance e/ ~ (is- Further, 75 = 2-kvs\Tsd\ 2 and 
7; = 2ttvi\Tdl\ 2 define the tunneling rates between su- 
perconductor and dots, and between dots and leads, re- 
spectively, with vs and V\ being the DOS per spin at 
the chemical potentials ^5 and (ii, respectively. We will 
show that the total effective tunneling rate from the su- 
perconductor to the leads is given by 75/7; due to the 
Andreev process. To specify the initial state \i) of the 
system we point out that since A > A/Lt, &_bT, the su- 
perconductor contains no quasiparticle initially. Also, 
due to the asymmetric barrier case, \Tdl\ > \Tsd\ ( or 
li > 7s) j an electron on the dot level e; leaves the dot 
to the leads much faster than new electrons can be pro- 
vided by the superconductor. In addition, a stationary 
occupation of the resonant levels e/ is given by the grand 
canonical distribution function oc exp(— A/i/ksT) which 
is exponentially small if A/i > fc^T. This implies pi ~ 
for initial states with occupied levels e;. Therefore we 
consider the initial state \i) — \0}s\0}d\hi)i, where |0)s 
is the quasiparticle vacuum for the superconductor, |0)d 
means that both dot levels e/ are unoccupied, and \(ii)i 
defines the occupation of the leads which are filled with 
electrons up to the chemical potential pi . We remark that 
in our regime of interest no Kondo effects appear which 
could destroy the spin entanglement, since our dots con- 
tain each an even number of electrons in the stationary 
limit. 



VI. STATIONARY CURRENT AND T-MATRIX 

The stationary current of two electrons passing from 
the superconductor via virtual dot states to the leads is 
given by 



2e^ 



f,i 



(6.1) 



where Wfi is the transition rate from the superconductor 
to the leads. We-,calculate this transition rate in a T- 
matrix approacf£3, 



W ft = 2Tr\(f\T(e t m 2 6(e f -e l 



(6.2) 



Here, T(e^) = Ht e . + l n _ H (si — Ho), is the on-shell trans- 
mission or T-matrix, with 77 being a small positive real 



VII. CURRENT DUE TO TUNNELING TO 
DIFFERENT LEADS 

We now calculate the current for the simultaneous tun- 
neling of two electrons into different leads. Since we as- 
sume that the spin is a good quantum number we can 
specify the final state for two electrons, one of them be- 
ing in lead 1 the other in lead 2, according to their total 
spin S. This spin can be either a singlet (in standard 
notation) |5) = (| U) - I 4T))/V2 with S = 0, or a 
triplet with S = 1. In the regime of interest (4.1), and 
since the total spin is conserved, [S 2 ,_ff] = 0, the sin- 
glet state of the initial Cooper pair will be conserved in 
the transport process and the final state must also be 
a singlet. That this is indeed true we can see explicitly 
if we allow for the possibility that the final state could 
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also be the S g = triplet \t ) = (| U) + | 4T))/\/2. The 
not-entangled triplets \t + ) = | ft) and |i_) = | J, J.) can 
be excluded right away since the tunnel Hamiltonian Ht 
conserves the spin-component a and an Andreev process 
involves tunneling of two electrons with different spin a. 
Therefore we consider two-particle final states of the form 



(i|(d 2 j.rfit ~ d 2^d H )T \i) 



in = (1/V2) (4 pT 4 



Ml ± a lp| a 2qT 



I*)- 



(7.1) 



The — and + signs belong to the singlet | S) and triplet 
| to), resp. Note that this singlet/triplet state is formed 
out of two electrons, one being in the p state in lead 1 
and with energy e p , while the other one is in the q state 
in lead 2 with energy e q . Thus, the two electrons are 
entangled in spin space while separated in orbital space, 
thereby providing a non-local EPR pair. The tunnel pro- 
cess to different leads appears in the following order. A 
Cooper pair breaks up due to tunneling of an electron 
with spin a to one of the dots (with empty level e/) from 
the point of the superconductor nearest to this dot. This 
is a virtual state with energy deficit > A. Since 
A > 7;, the second electron of the Cooper pair with 
spin — a tunnels to the other empty dot-level before the 
electron with spin a can escape to the lead. Therefore, 
both electrons tunnel almost simultaneously to the dots 
(within the uncertainty time h/A). Since we work at 
the resonance e; ~ /is = the energy denominators in 
( |6.3[ ) show divergences oc 1/t) indicating that tunneling 
between the dots and the leads is resonant and we have 
to treat tunneling to all orders in Hrji in (6.3). As a 
result r\ will be replaced by ji/2. Tunneling back to the 
superconductor is unlikely since \Tsd\ < \Tdl\- We can 
therefore write the transition amplitude between initial 
and final state as 

(f\ T o\i) = -j^(i|a 2q |ai pT T 4t4|H) 

x(*|(d ai di t ±d aT d u )T"|t) J (7.2) 

where To = T{ei — 0), and the partial T-matrices T and 
T are given by 



T — —H S d- -ttHsd 

it] - H ti] - H 



and 



T = H DL £ (—L-H DL ) 

^ V if) -Hq J 



n=0 



2n+l 



(7.3) 



(7.4) 



In (7.2) we used that the matrix element containing T is 
invariant under spin exchange T<-»J.. The part containing 
T describes the Andreev process of tunneling of two 
electrons with opposite spins from the SC to different 
dots 1,2, while the part containing T is the resonant dot 
«-> lead tunneling. 

We first consider the Andreev process. We insert a 
complete set of single-quasiparticle (virtual) states be- 
tween the two Hsd in ( |7.3| ) and find for the singlet final 
state 



ei + £2 



17] 



E 



£k 



cos(k-<5r), (7.5) 



where 8r = i"i — r 2 denotes the distance vector between 
the points on the superconductor from which electron 1 
and 2 tunnel into the dots. Note that the triplet contri- 
bution vanishes since itk^k — M_ kf- k for s-wave super- 
conductors. The sum over k in ( |7.5| ) can be calculated 
by linearizing the spectrum around the Fermi energy and 
we obtain 



E 



u k v k 



E> -.fr-*)-^^.-^, (7.6) 



where kp is the Fermi wavevector in the superconductor. 

VIII. DOMINANT RESONANT TUNNELING 
EVENTS 

We turn to the calculation of the matrix clement in 



(7.2) containing T where tunneling is treated to all or- 
ders in HrjL- We introduce the ket notation |12), and, 
for simplicity, suppress the spin index a. Here 1 stands 
for quantum numbers of the electron on dot 1 /lead 1 and 

similar for 2. For example, \pq) stands for £iJ po .a2q- CT K)> 
where p is from lead 1 and q from lead 2, or, correspond- 
ingly, \pD) stands for a\ pa d\ _ a \i), etc. We restrict our- 
selves to the resummation of the following dot <-> lead 
transitions \DD) -> \LD) -> \DD) or \DD) -> \DL) -> 
\DD). In this sequence, \DD) is the state with one elec- 
tron on dot 1 and the other one on dot 2, and \LD) de- 
notes a state where one electron is in lead 1 and the other 
one on dot 2. We thereby exclude tunneling sequences 
of the kind \DD) -► \LD) \LL) -> \LD) -> \DD) or 
\DD) -> \LD) \LL) -> \DL) -> \DD), where both 
electrons are virtually simultaneously in the leads as well 
as the creation of electron-hole-gair excitations out of the 
Fermi sea. We showed in Re f£Il that such contributions 
are suppressed in the regime ( f4.l| ) considered here by the 
small parameter ji/A/j.. It is quite clear that electron- 
hole pair excitations, which in principal could spoil the 
entanglement between the original partner electrons, are 
suppressed in our regime of interest. The following qual- 
itative argument is in order (for a detailed calculation, 
seetll). Suppose an electron initially on, say, dot 1, tun- 
nels to the lead 1. Instead of hopping back to the dot, 
thereby completing the sequence \DD) — » \LD) — > \DD), 
it disappears in a final state out in the lead (with energy 
£p ~ £ i)- Now another electron from the Fermi sea ,with 
energy e < — A/i, hops on the empty dot. The creation 
of such an electron-hole pair involves an energy deficit 
of at least A/x (for e\ =0). Since this is a virtual state 
it can only exist during the uncertainty time ~ 1/Afi. 
The relaxation of this energy deficit requires the tunnel- 
ing of this "Fermi sea" electron from the dot back to the 
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lead. Since the tunneling from the dot to the lead takes 
a time on the scale of I/71 this process is suppressed 
since A[i > 71. The dominant contributions are then 
resummed in the following sequence 

(pq\T'\DD) 

OO 

(pq\H DlLl \Dq)(Dq\ £ H DlLl ) 2n \Dq) 



71=0 



^{Dq\ T ^ Tro H D2L2 \DD) 

OO 

+{pq\H D2L2 \ P D)(pD\ £ (-Ijj-Hn^nD) 



71=0 



x(pD>|^tf Clil |-DD) 



00 

^ IT] - H 



I -D-D) 



(8.1) 



Since the sums for the transition |-DZ?) — * |-DD) via the 
sequences \DD) -> |LD) |DL>) and |DL>) -> \DL) -> 
1-D.D) are independent, we can write all summations in 



(3.1) as geometric series which allow for explicit resum- 
mations and we obtain 



(pq\T \DD) = 



T DL{ € i+ e 2-iri) 



*7i/2)(ea 



z 72 /2) 



(8.2) 



Thus, we see that the resummations in ( |8. ip cancel all 
diverge nces like the {e\ + e 2 — 277) denominator appear- 
ing in ( [7.5|) , and that, as expected, the resummation 
of divergent terms leads effectively to the replacement 
irj — ► iji/2 which makes the limit e; — * well-behaved. 
In ( |8.2| ) we neglected a small logarithmic correction to 
the (bare) resonant levels ej which is given as the real 
part of the self energy S ; = \T DL \ 2 £ k (w? - ej - £k) _1 
with |ReS;| ~ 7/ln(e c /A/i) < A/z and is therefore not 
important. The ener gy e e is the conduction band cut-off. 



Mak ing use of Eqs. (0,^2) and of Eqs. (7.5,7.6) and 
(3.2), we finally obtain for the current (denoted by I\) 
where each of-the two entangled electrons tunnels into a 
different leadO 



h = 



(ei + e 2 ) 2 +7 2 /4 



sm(kpdr) 
kp5r 



exp 



2Sr 



(8.3) 



where, 7 = 71 + 72. We note that Eq. ( |8.3| ) also holds 
for the case with 71 7^ 72 • The current becomes exponen- 
tially suppressed with increasing distance Sr between the 
tunneling points on the superconductor, the scale given 
by the Cooper pair coherence length £. This does not 
pose severe restrictions for conventional s-wave materi- 
als with £ typically being on the order of micrometers. 
More severe is the restriction that kpSr should not be too 
large compared to unity, especially if hp 1 of the super- 
conductor assumes a typical value on the order of a few 
Angstroms. Still, since the suppression in kpSr is only 
power-law like there is a sufficiently large regime on the 
nanometer scale for Sr where the current 1\ can assume 



a finite measurable value. The power law suppression of 
the current in 1/kpSr is very sensitive to the dimension 
of the SC and we suspect that the suppression will be 
softened by going over to lower dimensional supercon- 
ductors. We will address this issue in Section g. The 



current ( |8.3 ) has a two-particle Breit-Wigner resonance 
form which assumes it maximum value when e\ = — e 2 
(see also Fig. 3, and note that fxs = 0), 



h = 



4e7| 
7 



sin(kp5r) 
kf5r 



2Sr 



exp . (8.4) 



This resonance at ei = ~ti clearly shows that the current 
is a correlated two-particle effect (even apart from any 
spin correlation) as we should expect from the Andreev 
process involving the coherent tunneling of two electr ons . 
Together with the single-particle resonances in Eq. (|8.2| ) 
and by using energy conservation Ei — Sf = 0, which 
implies e p — — e q , we thus see that the current is carried 
by correlated pairs of electrons whose lead energies satisfy 



C-2 



< 



72- 



|e p - ei| ^5 7i and kc 

A particularly interesting case occurs when the ener- 
gies of the dots, ei and €2, are both tuned to zero, i.e. 
£1 = f-2 = us = 0. We stress that in this case the elec- 
tron in lead 1 and its spin-entangled partner in lead 2 
possess exactly the same orbital energy. We have shown 
previously!!-] that this degeneracy of orbital energies is 
a crucial requirement for noise measurements in which 
the singlets can be detected by an enhanced noise in the 
current (bunching) due to a symmetrical orbital wave- 
function of the singlet state, whereas uncorrelated elec- 
trons, or, more generally, electrons in a triplet state, lead 
to a suppression of noise (antibunching). Note that not 
all triplets are entangled states. Only the triplet with 
S z = is entangled. Measurement of noise enhancement 
is therefore a unique signature of entanglementJi-3. 

We remark again that the current I\ is carried by elec- 
trons which are entangled in spin space and spatially 
separated in orbital space. In other words, the station- 
ary current 1\ is a current of non-local spin-based EPR 
pairs. Finally, we note that due to the singlet character 
of the EPR pair we do not know whether the electron 
in, say, lead 1 carries an up or a down spin, this can 
be revealed only by a spin- measurement. Of course, any 
measurement of the spin of one (or both) electrons will 
immediately destroy the singlet state and thus the en- 
tanglement. Such a spin measurement (spin read-out) 
can be performed e.g. by making use of the spin fil- 
tering effect of quantum dotsca. The singlet state will 
also be destroyed by spin-dependent scattering (but .apt 
by Coulomb exchange interaction in the Fermi sealij). 
However, it is known experimentally that electron spins 
in a semiconductor environment show unusually long 
dephasing times approaching microseconds and can be 
transported phase coherently over distances exceeding 
100/i/TUjiD. This distance is sufficiently long for experi- 
ments performed typically, on the length scale of quantum 
confined nanostructuresc-3. 
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IX. CURRENT DUE TO TUNNELING VIA THE 
SAME DOT 

In this section we calculate the current for tunneling of 
two electrons via the same dot and subsequently into the 
same lead. We show that such processes are suppressed 
by a factor ("fi/U) 2 and/or (7;/ A) 2 compared to the pro- 
cess discussed in the preceding section. But, in contrast 
to the previous case, we do not get a suppression result- 
ing from the spatial separation of the Cooper pair on the 
superconductor, since here the two electrons tunnel from 
the same point either from ri or r2 (see Fig. 2). As be- 
fore, a tunnel process starts by breaking up a Cooper pair 
followed by an Andreev process with two possible tunnel- 
ing sequences, see Fig. 4, panel a). In a first step, one 
electron tunnels from the superconductor to, say, dot 1, 
and in a second step the second electron also tunnels to 
dot 1. Now two electrons are simultaneously on the same 
dot which costs additional Coulomb repulsion energy U , 
thus this virtual state is suppressed by l/U. Finally, the 
two electrons leave dot 1 and tunnel into lead 1. There 
is an alternative competing process, see Fig. 4, panel b), 
which avoids the double occupancy. Here, one electron 
tunnels to, say, dot 1, and then the same electron tunnels 
further into lead 1, leaving an excitation on the super- 
conductor which costs additional gap energy A (instead 
of U) , before finally the second electron tunnels from the 
superconductor via dot 1 into lead 1. 

We first concentrate on the tunneling process b), and 
note that the leading contribution comes from the pro- 
cesses where both electrons have left the superconductor 
so that the system has no energy deficit anymore. We 
still have to resum the tunnel processes from the dot to 
the lead to all orders in the tunnel Hamiltonian Hol- 
In what follows we suppress the label I = 1,2 since the 
setup is assumed to be symmetric and tunneling into ei- 
ther lead 1 or lead 2 gives the same result. The transition 
amplitude (f\To\i) including only leading terms is 

(/|2b|i) = 

OO 1 

]T(/|H DL J2(-- TrH DL rWa) 

L — ' z — ' ir\ — Hq 



:(Dp"a\ 



71=0 
1 



-H, 



1 



sv 



-H, 



1 



DL " 



in - H Q ir) - Ho in - Ho 



-H SD \i), (9.1) 



± vlvt)!*)' with 



where again |/) = (1/V2) (a^a^ 
± denoting the triplet (+) and singlet (— ), resp., and 
the intermediate state \Dp"a) = d}_ a a L, There are 



some remarks in order regarding Eq. (9.1). The electron 
which tunnels to the state |p"er) has not to be resummed 
further since this would lead either to a double occupancy 
of the dot which is suppressed by 1//7, or to the state with 
two electrons simultaneously in the lead with a virtual 
summation over the state p". But we already mentioned 
above that the latter process is suppressed by ii/A/j,. 
The first factor in ([O]) describes therefore the multiple 





FIG. 4. Two competing virtual processes are shown when 
the two electrons tunnel via the same dot into the same lead: 
a) Andreev process leading to a double occupancy of the dot 
with virtual energy suppression l/U , and b) the process which 
differs by the sequence of tunneling, leading to an additional 
virtual energy suppression 1/A instead of l/U. 



dot <-> lead tunneling of the electron with spin — a which 
resides on the dot in the intermediate state \Dp"a) = 
d- a a p , lcr \i) and which eventually tunnels in its final state 
in the lead. Again, this amplitude can be resummed 
explicitly with the result for a =f, 



1 ill - H 
Tdl ei + e F 



IT) 



\/2 ei + e 
and for a =| we get 

1 



ili/2 



( V' P T S] 



p"p' j 



(9.2) 



(WdlYx 



71 = ' 



Ho 

Tdl e; + e P 



H DL ) 2n \Dp" [) 



V2 ei+ e p » - iji/2 



(WtW> (9-3) 



In (9.2) and ( |9.3| ) the upper sign belongs again to the 
triplet and the lower sign to the singlet. For the am- 
plitude containing the superconductor-dot transitions in 
(|9.lD we obtain 
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(Dp"a 



1 1 1 

Hsd- ttHdl- —HsdM 



ir] - H ir) - H ir] - H Q 
sT DL Tl D v s 



A(e, 



-ii])' 



(9.4) 



wher e s = +1 (-1) for a =| {{). Combining the results 
(|]2])-(Q we obtain for the amplitude 



(/|3b|i) 



2 3 / 2 v s (T SD T DL ) 2 (ei-i 7l /2) 



A(e, 



*Ti/2)(Q 



V 



i-yi/2) 



(9.5) 



for the singlet final state, whereas we get again zero for 
the triplet. 

Next we consider the process where the tunneling in- 
volves a double occupancy of the dot, see panel a) in Fig. 
4. In this case the transition amplitude is 

(/|To|i) = 

oo 1 

p"cr " '' ' 



:(Dp"a\ 



n=0 
1 



-H DL 



1 



H S d- 



1 



ir\ -Ho ir]- Ho ii] - H 



-HsdH) 



(9.6) 



Repeating a similar calc ulat ion as before we find that the 
amplitude is given by ( |9.5| ) but with A being replaced 
by U/ir, and again, (f\T \i) is only nonzero for the sin- 
glet fina l state. We note that the two amplitudes (9.5) 
and (| 



^ have the same initial and same final states and 
therefore have to be added coherently to obtain the total 
current due to processes a) and b). Then, using Eq. (6.1) 
we find for the total current I?, in case of tunneling of two 
electrons into the same leadtH 



e7s7 
£ 2 : 



1 

I 



1 

7rA 



1 

u 



(9.7) 



We see that the effect of the quantum dots shows up 
in the suppression factor (j/£) 2 for tunneling into the 
same lead. We remark that in contrast to the previous 
case (tunneling into different leads) the current does not 
have a resonant behavior since the virtual dot states are 
no longer at resonance due the energy costs U or A in 
the tunneling process. We now compare I\ given in (8.4) 
with Ii by forming the ratio of the two currents 



h 
h ' 

From this 
dominatin; 
We would 
sion of h 
rent I\) is 

e; - Ms = 
of the two 
achieved. 



AS 2 



sm(kpSr) 



kpSr 



exp 



2Sr 



(9.8) 



ratio we see that the desired regime with I\ 
1 1% is obtained when S/'y > kpSr, and Sr < £. 
like to emphasize that the relative suppres- 
(as well as the absolute value of the cur- 
maximized by working around the resonances 
= f£j and, in addition, the desired injection 
electrons at the same orbital energy is then 



X. EFFICIENCY AND DISCUSSION 

The current I\ and therefore the ratio ( |9.S| ) suffers an 
exponential suppression on the scale of £ if the tunneling 
of the two (coherent) electrons takes place from different 
points i"i and Y2 of the superconductor. For conventional 
s-wave superconductors the coherence length £ is typ- 
ically on the order of micrometers and therefore poses 
not severe restrictions. So in the interesting regime the 
suppression of the Andreev amplitude is only polynomial 
oc 1/kpSr. It was shownE3 that a superconductor on top 
of a two-dimensional electron gas (2DEG) can induce su- 
perconductivity (by the proximity effect) in the 2DEG 
with a finite order parameter. The 2DEG then becomes 
a two-dimensional (2D) superconductor. One could then 
desire to implement the two quantum dots in the 2DEG 
directly. More recently, it was suggested that supercon- 
ductivity should also be present in ropes of single- walled 
carbon nanotubesEj which are strictly one-dimensi ona l 
(ID) systems. It is therefore interesting to calculate (7.6) 
also in 2D and ID. In the case of a 2D superconductor 
we evaluate (|7.6|) in leading order in Sr / 7r£ and find 



-g— cos ( k • 

k(2£>) ^ k 



% [Jo(k F Sr)+2j2 



(10.1) 



v=l 



The right-hand side of (10.1) can be approximated by 
(it/2)vs Jo{kF5r){\ — (2/ir) ln2) in the limit of large kpSr. 
For large kpSr, the behavior of the zeroth-order Bessel- 
function is Jo(kpSr) ~ y/2/7rkpSr cos(kp Sr — (tt/4)). So 
the amplitude decays asymptotically only oc 1/ \/kpSr or 
the current I\ by a factor cx 1/kpSr, respectively. In the 
case of ID we obtain 



E 

k(ir>; 



cos (k ■ Sr) 



% s cos{k F Sr) e-W*® , 



(10.2) 



where there are only oscillations and no decay of the 
Andreev amplitude (for 8r/x£ < 1). We see that the 
suppression due to the finite separation of the tunnel- 
ing points on the superconductor can be reduced con- 
siderably (or even excluded completely) by going over to 
lower-dimensional superconductors. By taking into ac- 
count the dependence on the dimension of the supercon- 
ductor we can relax the condition for the entangler to be 
efficient to 



(£/j) 2 > (k F Sr) d - 



(10.3) 



where d is the dimension of the superconductor. 
We note that the coherent injection of the two spin- 
entangled electrons by an Andreev process via the dots 
into the leads allows for a time resolved detection of indi- 
vidual Cooper pairs in the leads since the delay time be- 
tween the two partner electrons of a Cooper pair is given 
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by 1/A whereas the separation in time of subsequent 
Cooper pairs is given approximately by 2e/l\ ~ 111 Is- 
Since A > 75 and, in addition, 7/ > 7s the time delay 
between the two partners of the original Cooper pair is 
much shorter than the time difference between different 
Cooper pairs. 



XI. AHARONOV-BOHM OSCILLATIONS 

In this section we show that the different tunneling 
paths of the two electrons from the superconductor to the 
leads can be detected via the flux-dependent Aharonov 
Bohm oscillations in the current flowing through a closed 
loop (see Fig. 5). We show that due to the possibil- 
ity that two electrons can tunnel either via different 
dots into different leads (non-local process) or via the 
same dot into the same lead (local process), the cur- 
rent as a function of magnetic flux ip penetrating the 
loop contains h/e and h/2e oscillation periods. To be 
concrete, we consider a setup where the two leads 1 
and 2 are connected such that they form an Aharonov- 
Bohm loop, (see Fig. 5), where the electrons are in- 
jected from the left via the superconductor, traversing 
the upper (lead 1) and lower (lead 2) arm of the loop be- 
fore they rejoin to interfere and then exit into the same 
lead, where the current is then measured as a function 
of varying magnetic flux <f>. In the presence of a mag- 



d 2 l 2 



sc 




FIG. 5. The setup where the two outgoing leads 1 (LI) and 
2 (L2) are connected to a common lead so that the tunneling 
path of the electrons can form a loop. By applying a magnetic 
flux <f> the current shows Aharonov-Bohm oscillations with 
periods h/e and h/2e which can be used to identify different 
tunneling paths of the two electrons. 

netic flux, each tunneling amplitude obtains a phase fac- 



tor, T DlLl -> T DlLl e^' 2 *\ and T< 



d 2 l 2 



) 

where (f>o = h/e is the single-electron flux quantum. For 
simplicity of the discussion we assume that the entire 
phase is acquired when the electron hops from the dot 
into the leads, so that the process dot-lead-dot gives 
basically the full Aharonov-Bohm phase factor e ±Ic ^ 
of the loop (and only a negligible amount of phase is 
picked up along the path from the superconductor to the 
dots). We stress that there is no loss of generality in this 
assumption. The transition amplitude from the initial 
state to the final state has now the following structure 



(f\T \i) ~ T DlLl T D2L2 + T^jJ+lt* + T\ 
Here, the first term comes from the process via differ- 
ent leads (see ( |8.2| )), where no Aharonov-Bohm phase 
is picked up. The Aharonov-Bohm phase appears in 
the remaining two terms, which come from processes 
via the sa me leads, either via lead 1 or lead 2 (see 
( |9.5| ) and ( |9.6|) ). The total current / is now obtained 
from I (/|To|i) | 2 together with a summation over the final 
states, giving I = I%+ 12 + Iab, and theJlux-dependent 
Aharonov-Bohm current Iab is given byEZ 



Iab = y/8hI 2 F(ei)cos 



■o)+J 2 cos(2^o), (11.1) 

£l =, (11.2) 

V^ + (7l/2) 2 ' 



where, for simplicity, we have assumed that e\ = e 2 = (4, 
and 71 = 72 = 7l- Here, the first term (different leads) is 
periodic in fig like for single-electron Aharonov-Bohm in- 
terference effects, while the second one (same leads) is pe- 
riodic in half the flux quantum <fro/2, describing thus the 
interference of two coherent electrons travelling the up- 
per or the lower arm of the loop (similar single- and two- 
particle Aharonov-Bohm effects occur iy-the Josephson 
current through an Aharonov-Bohm loopEjo) . It is clear 



from (11.1) that the h/e oscillation comes from the inter- 
ference between a contribution where the two electrons 
travel through different arms with contributions where 
the two electrons travel through the same arm. Both 
Aharonov-Bohm oscillations with period h/e, and h/2e, 
vanish with decreasing I2, i.e. with increasing on-site re- 
pulsion U and/or gap A. However, their relative weight 
is given by \J T\ /I 2 , implying that the h/2e oscillations 
vanish faster than the h/e ones. This behavior is quite 
remarkable since it opens up the possibility to tune down 
the unwanted leakage process ~ I2 cos (20/<fo) where two 
electrons proceed via the same dot/lead by increasing U 
with a gate voltage applied to the dots. The dominant 
current contribution with period h/e comes then from the 
desired entangled electro ns pro ceeding via different leads. 
On the other hand, if */ h /I2 < 1, which could become 
the case e.g. for UfSt > £/j, we are left with h/2e oscil- 
lations only. Besides the fact that the Aharonov-Bohm 
oscillations are interesting in its own right, the Aharonov- 
Bohm oscillations further provide an experimental probe 
of the non-locality of the two spin-entangled electrons. 
Note that dephasing processes which affect the orbital 
part suppress Iab- Still, the flux-independent current 
I\ + 12 can remain finite and contain electrons which are 
entangled in spin-space, provided that there is only neg- 
ligible spin-orbit coupling so that the spin is still a good 
quantum number. 

We would like to mention another important feature of 
the Aharonov-Bohm effect under discussion, namely the 
relative phase shift between the amplitudes of tunneling 
to the same lead and to different leads, resulting in the 
additional prefactor F(ci) in the first term of the right- 
hand side of Eq. (11.1). This phase shift is due to the 
fact that there is a two-particle resonance in the ampli- 
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tude ( |3.2j ) while th ere is onl y a single-particle resonance 
in the amplitudes (9.5) and (9J3) (we recall that the sec- 
ond resonance is suppressed by the Coulomb blockade 
effect). Thus, when the chemical potential /is of the 
superco nduc tor crosses the resonance, |e;| < 7^, the am- 
plitude (|8.2|) acquires an extra phase factor e l ^ r , where 
4> r = arg[l/(e; — vyi,/2)]. Then the interference of the 
two amplitudes leads to the prefactor F ( e;) = cos<p r in 
the first term on the right-hand side of flll.l| ). In par- 
ticular, exactly at the middle of the resonance, e; = 0, 
the phase shift is <f> r — tt/2, and thus the h/e oscillations 
vanish, since F(0) = cos(7r/2) = 0. Note however, that 
although F = ±1 away from the resonance (|e;| > jl) 
the h/e oscillations vanish again, now because the current 
ii ~ e7|7L/ef vanishes. Thus the optimal regime for the 
observation of the Aharonov-Bohm effect is |e;| ~ 7^. 

Finally, the preceding discussion shows that even if 
the spins of two electrons are entangled their associ- 
ated charge current does not reveal this spin-correlation 
in a simple Aharonov-Bohm interference experiments. 
Only if we consider the current-current correlations 
(noise) in a beam splitter setup, can we detect also this 
spin-correlation in the transport current via its charge 
properties^. 



XII. ANDREEV ENTANGLER WITH 
LUTTINGER LIQUID LEADS 

We discuss a further implementation of an entangler 
for electron spins. We replace the quantum dots and the 
noninteracting leads by interacting one-dimensional wires 
described as Luttinger liquids (LL) (see Fig. 6). The low 
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FIG. 6. The entangler setup with Luttinger liquid leads: 
An s-wave superconductor, SC, with chemical potential fis 
is weakly coupled to two strongly interacting bulk Luttinger 
liquids 1,2 held at the same chemical potential (Xi. The two 
spin-entangled electrons of a Cooper pair can coherently tun- 
nel with amplitude to by means of an Andreev process from 
the superconductor to the leads. Two competing transport 
channels are considered. Two electrons can tunnel to differ- 
ent leads from points ri and r2 (with distance Sr) of the SC 
or from the same point ri or V2 into the same lead. The 
interaction between the two leads is assumed to be negligible. 

energy excitations of these LL are collective charge and 
spin density oscillations rather than quasiparticles which 
resemble free electrons. As a consequence tunneling into 



a LL is strongly suppressed due to interaction in the LL- 
leads. We show that the interaction can be used to sepa- 
rate two electrons, originating from an Andreev process, 
so that they preferably tunnel into different leads rather 
than into the same. Again we take into account a fi- 
nite tunneling distance on the superconductor when the 
electrons tunnel to different leads. 



XIII. MODEL 

We consider an s-wave superconductor as described in 
Section ^ with chemical potential /is which is tunnel- 
coupled to two (spatially separated) bulk LL-leads, both 
held at the same chemical potential /i/ so that a finite bias 
fi = pbs — fj,i is applied. The Hamiltonian of the whole 
system is H = H + H T , with H = H s + J2 n =i,2 H Ln 
describing the isolated superconductor and LL-leads 1,2, 
resp.. Tunneling between the superconductor and the 
leads is governed by the tunnel Hamiltonian Ht, which 
will be described below. 

The two leads 1, 2 are assumed to be infinitely extended 
and interacting one-dimensional systems described by 
conventional LL-theory. The LL-Hamiltonian for the 
low energy excitations of lead n — 1,2 is written in a 
bosonized form as (neglecting backscattering)E3 



+L/2 



v =P' a -L/2 



Hhn — HlN n — 



^^nL + ^-(^) 2 ). (18.1) 



where N n is the number operator for lead n = 1,2, 
and the fields Tl n (x) and 4> n {x) satisfy bosonic commu- 
tation relations \^ nv (x),U mfl (x')} = i5 nm 5 upi 5{x - x 1 ). 



The Hamiltonian (13.1) describes long- wavelength charge 
(y = p) and spin [v = a) density oscillations in the 
LL propagating with velocities u p and u a , resp. The 
velocities u v and the stiffness parameters K v depend 
on the interaction between the electrons in the LL. In 
the limit of vanishing backscattering, u a — vp and 
K a = 1 and the LL is described by only two param- 
eters K p < 1, and u p = Vp/K p , with vp being the 
Fermi velocity. We decompose the field operator describ- 
ing electrons with spin s into a right and left moving 
part, ip ns {x) = e lk F*ip ns+ (x) + e^ k F*ip ns _( x ). The right 
(left) moving field operator ip ns+ (x) (ip ns - (aj)Iis-then ex- 
pressed as an exponential of bosonic fields astjO 

i>ns±{x) = lim -±2gl exp ( ± — ((j> np {x) 

V27TQ! I- V2 ^ 

+ S(f) na (x) T (&np(x) + Sdncr{x))j j 

(13.2) 

where d x 6 nv (x) = TrTl nu (x). The operators r)± ins are 
needed to ensure the correct fermionic anticommutation 
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relations. In (13.2) and hereafter we adopt the conven- 
tion that s = +1 for s =|, and s = —1 for s =1, if s has 
not the meaning of an index. 

Transfer of electrons from the SC to the LL-leads 
is described by the tunnel Hamiltonian Ht — 
*o J2 na V4 s *s( r n) + H.c. The operator * s (r n ) anni- 
hilates an electron with spin s at the point r„ on the 
superconductor nearest to the LL-lead n, and ij}\ a cre- 
ates it again with amplitude to at point x n in LL n (see 
Fig. 6). We assume that the tunneling amplitude t does 
not depend on spin and also is the same for both leads. 



XIV. CALCULATION OF THE CURRENT: 
TWO COMPETING CHANNELS 



suppressed if fi < A with the result, in leading order in 
f/A, 



b=±i 



(14.3) 



The constant Af, in (14.3) is of order one for not too 
strong interactions, but is decreasing with increasing in- 
teractions in the LL-leads, and is given by 



A h = - 



r(2 7p + 2) 



r(7 P h + 1) 



2 r(2 7 p b + 2 lp + 1) i \ r(ig + i)r(i^) / 



We now calculate the current for tunneling of two spin- 
entangled electrons into different and into the same leads. 
Again, we use a T-matrix approach introduced in Sec- 
tion VI and calculate the current in lowest order in the 



tunneling Hamiltonian Ht, which describes an Andreev 
process. For the current I\ for tunneling of two electrons 
into different LL-leads we obtain in leading order in fi/A 
and at zero temperature 



h = 



T(2 7p + 2) u p 



2A/i 



27p 



(14.1) 



where T is the Gamma function and A is a short dis- 
tance cut-off on the order of the lattice spacing in the 
LL. The interaction suppresses the current considerably 
and the bias dependence has its characteristic non-linear 
form Ii oc /i 27 p +1 with an interaction dependent expo- 
nent 7 P = (K p + Kp 1 )/-!- 1/2 > 0. The parameter 7 P is 
the exponent for tunneling into the bulk of a single LL, 
i.e. p(e) ~ |e| 7p , where p(e) is the single particle DOSE3. 
The noninteracting limit I® is given as 



1° = Trej 2 fj, 



sin(/cp(5r) 



kpSr 



25r 



exp , (14.2) 



with 7 = Airvgvi |io| 2 being the probability per spin to 
tunnel from the SC to the LL-leads and v$ and v\ are 
the energy DOS per spin for the superconductor and the 
LL-leads at the chemical potentials /is and \ii resp. . We 
again observe the same dependence on the tunneling dis- 
tance Sr as before in the Andreev Entangler (see flS.3p). 
The results for the lower dimensional superconductors 
are as before and given in Section [x|. 
Now we compare this result with the process when the 
two electrons tunnel into the same lead (with Sr = 0). 
Note that in this case the correlation in time between sub- 
sequent tunneling events of the two electrons of a Cooper 
pair becomes extremely important due to the interaction 
in the LL-lead. We find after some calculation that the 
current I2 for tunneling into the same lead (1 or 2) is 



(14.4) 



We remark that in (14.3) the current I\ is to be taken 
at Sr = 0. The noninteracting limit I2 = h = I± is 
rediscovered by putting j p — 7/3 & = 0, and u p = Vf- 
The result for I2 shows that the unwanted injection of 
two electrons into the same lead is suppressed compared 
to Ii by a factor of (2^/A) 27p + , where -f p+ = j p , if 
both electrons are injected into the same branch (left 
or right movers), or by (2/i/A) 27p - if the two elec- 
trons travel in different directions. Since it holds that 
7 P - = 7p+ + (1 — K p )/2 > 7 P +, it is more favorable 
that the two electrons travel in the same direction than 
that in opposite directions. The suppression of the cur- 
rent I2 by 1/A shows very nicely the two-particle cor- 
relation effect in the LL, when the electrons tunnel into 
the same lead. The larger A, the shorter is the delay 
time between the arrivals of the two partner electrons of 
a Cooper pair, and, in turn, the more the second electron 
tunneling into the same lead will feel the existence of the 
first one which is already present in the LL. By increas- 
ing the bias \x the electrons can tunnel faster through 
the barrier since there are more channels available into 
which the electron can tunnel, and therefore the effect 
of A is less pronounced. Also note that this correlation 
effect disappears when interactions are absent (i.e. when 
7 p = 7pb = 0) in the LL. Actual experimental systems 
which show LL-behavior are e.g. metallic carbon nan- 
otubes with similar exponents as derived hereEaO. The 
exponents for tunneling into the bulk of a LL and for 
tunneling into the end of a Lb-aee. in general, different, 
and it is predicted by theoryE3~ll3 and consistent with 



experimentErc3, that ^buik < jend_ therefore expect 
that if the electrons tunnel into the end of the LL, the 
suppression is even more pronounced, but to find the 
correct exponents would require a careful recalculation 
in that geometry. 

Finally, the entangler setup is efficient if approximately 



27„+ 



> (k F Sr) c 



(14.5) 
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where d is the dimension of the superconductor, kp is 
the Fermi wavevector of the superconductor, and it is 
assumed that the coherence length £ of the SC is large 
compared to Sr (see Section |x|). 



XV. DECAY OF THE ELECTRON-SINGLET DUE 
TO LL-INTERACTIONS 

Here we consider the decay of a nonlocal singlet where 
one electron is in lead 1 and the other in lead 2 due to the 
interactions in the LL- leads. A quantity which accounts 
for this consideration is the correlation function 



P(v,t) = \(S(T,t)\S(0,0))\ 



(15.1) 



Here, P(r, t) is the probability that a singlet state in- 
jected at t — and at point r = (x\, x 2 ) = is found at 
some later time t at point r, and it is therefore a mea- 
sure of how much of the initial singlet state remains after 
switching on the interaction during the time interval t. 
The singlet state created on top of the LL-groundstates 
is 



|S(r,f)}=vM$ T (zi,t)^J ; (a:2,f) 
-^ n {x 1 ,t)^{x 2 ,t))\0), 



(15.2) 



where the extra normalization factor \Z2na for the singlet 
is introduced to guarantee J dr P(r, t) = 1 in the nonin- 
teracting limit. The singlet-singlet correlation function 
factorizes into a product of two single-particle correla- 
tion functions due to negligible interactions between the 
two leads 1,2 and we obtain for P(r,t), retaining only 
nonoscillatory terms, 



(15.3) 



with a time decaying weight factor of the 5-function 



f ® = n 



A 2 



A 2 + (v F - u v )H 2 



A 4 



\ (A 2 + (vFt) 2 - (u v t) 2 ) 2 + {2Ku v t) 2 ) 

(15.4) 

Without interaction we have F(t) = 1 which means that 
there is no decay of the singlet state. As interactions are 
turned on we see that for times t > A/u p the singlet state 
starts to decay in time. For long times t and for u a = vf, 
K a = 1 and u 

decay is F(t) 

very strong interactions in the LL leads, i.e. n p 
smaller than one, b ecom es F(t) ~ [A/it p i] 27p+1 . This 
result together with (15.3) shows that charge and spin of 
an electron propagate with velocity vp, whereas charge 
(spin) excitations of the LL propagate with u p (it<x) as 
we confirm explicitly in the next Section. 



p — vf/K p , the asymptotic behavior of the 

n k K J ,n A \^ V p [j} 2lp+ \ which for 
Up{X-K p ) L «2(1-^)J L t J 7 

K„ much 



XVI. PROPAGATION OF CHARGE AND SPIN 

Although the singlet state gets destroyed due to in- 
teraction the spin information can still be transported 
through the wires by the spin density excitations as we 
will show now. The normal ordered charge density p n {x) 
of lead n is p n (x) = J2 sr '■ ^nar WVw 0*0 where we 
kept only the slow spatial variations of the density op- 
erator, and r = ± denotes the branch of left- and right- 
movers. Likewise, for the normal ordered spin density 
we have <J z n (x) = J2sr s '■ ^LrW^W^) : • We now 
consider a state \^} = tp^ sr (x n )\0} where we inject an 
electron with spin s into branch r of lead n, and at time 
t = 0, on top of the LL groundstate. We calculate the 
time dependent charge and spin density fluctuations for 
this state and find for the charge density fluctuations 

(2?ra) (0\ip nsr (x n )p n (x' n , t)%l)\ sr {x n )\Q) 

= ^C 1 + rK p) S ( x 'n ~x n - u p t) 

+±(l-rK p )5{x' n -x n + u p t), (16.1) 

and for the spin density fluctuations 

(2na) (0| (a:„K « , t)^ nsr {x n )\0) 
= -(1 + rK a )8 (x' n -x n - u a t) 

+|(1 - rK a )5 (x' n - x n + U<7 t) . (16.2) 

We see that in contrast to the singlet, the charge and 
spin density fluctuations in the LL created by the in- 
jected electron do not decay and show a pulsed shape 
with no dispersion in time. This is due to the linear en- 
ergy dispersion relation of the LL-model. In carbon nan- 
otubes such a highly linear dispersion relation is indeed 
realized and therefore carbon nanotubes should be well 
suited for spi n tra nspor t. An other interesting effect that 
shows up in ( 16.1 ) and ( |16.2[ ) is the different velocities of 
spin and charge, which is known as spin-charge separa- 



tion. It would be interesting to test Bell incqualitiesE3 via 
spin-spin correlation measurements between the two LL- 
leads and see if the initial entanglement of the spin singlet 
is still observable in the spin density-fluctuations. Al- 
though detection of single spins with magnitudes on the 
order of electron spins has still not been achieved, mag- 
netic resonance force micrpacopy (MFRM) seems to be 
very promising in doing soc3. Another scenario is to use 
the LL just as an intermediate medium which is needed to 
first separate the two electrons of a Cooper pair and then 
to take them (in general other electrons) out again into 
two (spatially separated) Fermi liquid leads where the 
(possibly reduced) spin entanglement could be measured 
via the current noise in a beamsplitter experiments. In 
this conte xt w e note that the decay of the singlet state 
given by ( |15.3|) sets in almost immediately after the in- 
jection into the LLs (the time scale is approximately the 
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inverse of the Fermi energy) but at least at zero temper- 
ature, the suppression is only polynomial in time which 
suggests that some fraction of the singlet state can still 
be recovered. 



XVII. CONCLUSION 

We discussed setups for creating spin-entanglement 
in quantum confined nanostructures like semiconduc- 
tor quantum dots and mesoscopic wires based on su- 
perconductors which provide a natural source of spin- 
entanglement in form of Cooper pairs. We showed 
that a superconductor can induce tunable nonlocal spin- 
entanglement between two electron spins which reside on 
different quantum dots without having a direct tunnel 
coupling between the dots. We then described an en- 
tangler device that creates mobile spin-entangled elec- 
trons via an Andreev process into different dots which 
are tunnel-coupled to leads. The unwanted process of 
both electrons tunneling into the same leads can be sup- 
pressed by increasing the Coulomb repulsion on the quan- 
tum dots. We have shown that there exists a regime of 
experimental relevance where the current of entangled 
electrons shows a resonance and assumes a finite value 
with both partners of the singlet being in different leads 
but having the same orbital energy. This entangler thus 
satisfies the requirements needed to detect spin entan- 
glement via transport and noise measurements. Further, 
we discussed the flux-dependent oscillations of the cur- 
rent in an Aharonov-Bohm loop which could serve as an 
experimental means to detect the nonlocality of the two 
correlated electrons. Finally, we discussed a novel setup 
consisting of two interacting Luttinger liquid leads (such 
as nanotubes) weakly coupled to a superconductor to sep- 
arate the two spin-entangled electrons. We found that 
the coherent subsequent tunneling of two electrons into 
the same LL is suppressed if the applied bias is smaller 
than the superconducting gap in a characteristic power- 
law manner. 
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